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Abstract 

Supersymmetric models often predict a lightest superpartner (LSP) which is 
electrically charged and stable on the timescales of collider experiments. If such 
a particle were to be observed experimentally, is it possible to determine whether 
or not it is stable on cosmological timescales? Charged, stable particles are usu- 
ally considered to be excluded by cosmological arguments coupled with terrestrial 
searches for anomalously heavy water molecules. But when the cosmology is sig- 
nificantly altered, as can happen in models with large extra dimensions, these ar- 
guments are in turn significantly weakened. In this paper we suggest an alternate 
way to use searches for superheavy water to constrain the lifetimes of long-lived, 
charged particles, independent of most cosmological assumptions. By considering 
SUSY production by cosmic rays in the upper atmosphere, we are able to use cur- 
rent bounds on superheavy water to constrain the mass scale of squarks and gluinos 
to be greater than about 230 GeV, assuming a stable, charged LSP. This bound 
can be increased, but only by significantly increasing the size of the initial water 
sample tested. 



1 Introduction 



One of the primary goals of the current generation of colliders is the discovery of the 
mechanism by which the mass scale of the weak interactions is generated and stabilized. 
Whether that mechanism is supersymmetry (SUSY), technicolor, extra dimensions or 
something not yet considered, it will generically require a number of new particle states. 
One or more of these particles are often stable, at least on experimental time scales, 
thanks to discrete symmetries in the model or suppressed couplings. The archetypes for 
such particles are the "LSP" (lightest SUSY particle) in the most SUSY models, and the 
NLSP (next-to-lightest SUSY particle) in gauge-mediated models. 

If a new particle produced at a collider lives long enough to escape the detector, 
there is no simple way to measure its lifetime. Of utmost importance will be determining 
whether the lifetime is indeed finite or if the particle is stable on cosmological timescales. 
It has become standard lore that any stable remnant of new weak-scale physics must be, at 
most, weakly interacting. The possibility that the stable particle is electrically charged or 
strongly interacting has long been dismissed based on a combination of astrophysical and 
cosmological arguments. If such a particle, often called a CHAMP M, exists in abundance 
in our galaxy, then a steady flux of CHAMPs would have been captured by the earth over 
the course of our planet's lifetime. Such stable, charged particles would form hydrogen- 
like atoms which could in turn form superheavy water molecules. These superheavy 
atoms and molecules would collect in our oceans in significant abundances. Searches for 
superheavy isotopes of water in ocean water have all yielded negative results [0, [|. These 
searches have effectively ruled out charged dark matter at the weak scale. 

This result is obviously dependent on the flux of CHAMPs onto the earth. There 
are two theoretical avenues for determining this flux. The first assumes that these relics 
represent the majority component of the galactic dark matter halo; then their local 
abundance can be estimated directly from astrophysical data. From this we easily deduce 
the rate at which the CHAMPs would be captured and accumulate in our oceans. 

The second approach is to use standard cosmological assumptions and techniques to 
calculate the relic abundance of the CHAMPs from first principles. Then using calcula- 
tions of dark matter clumping, a local density can be estimated. These calculations also 
predict far greater abundances of superheavy water than are observed experimentally, 
even when the abundance is too small to account for the dark matter. 

There has been recent interest in whether charged stable relics are ruled out in more 
cosmologies. Motivated in part by theories with large extra dimensions studies have 
been done to calculate relic abundances for the case of extremely low reheat temperatures 
(as low as an MeV). The surprising result of these studies is that particles with TeV 
masses can still be created during the reheat process even if the final temperature is 
around the MeV scale [Q. When applied to electrically charged relics, a bound of roughly 
1 TeV on the relic's mass can be deduced using experimental search bounds |J. 

But can we bound, exclude or search for stable, charged relics without recourse to 
cosmology? In some models, physics near the TeV scale is so non-canonical that it is 
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unclear whether existing cosmological bounds apply at all (see for example the model 
of Ref. 0). In this paper, we will approach the same problem but from an original 
point of view. We will consider charged relics created by cosmic rays in the earth's 
upper atmosphere. In order to be specific, we will consider SUSY as our template model; 
however extension of these bounds to any other model of weak-scale physics with a stable 
charged relic would be trivial. 

We will not place a bound on the mass of the SUSY relic itself. Instead we will 
place a bound on the SUSY scale, or more specifically, the masses of the squarks and 
gluinos. Direct production of the charged relic in cosmic ray collisions may be very small 
(for example, the relic could be a stau with only small electromagnetic cross-sections). 
However, production of any SUSY state would, after a cascade of decays, result in one or 
more of the charged LSPs. Thus the production cross-section will generally be dominated 
by production of squarks and gluinos, not the LSP itself. 

None of these results depend in any way on new models which attempt to explain 
cosmic rays at the very highest energies. These results are generated using the usual 
spectrum of incident protons and heavy nuclei. Our only assumption is that this spectrum 
has remained essentially unchanged over the last 3 billion years of the earth's history. 

2 Cosmic Rays as Particle Colliders 

While the energy frontier in particle accelerators is slowing moving higher and higher, 
collisions with center-of-mass energies in the TeV range have been common throughout 
the history of the earth in the form of cosmic rays scattering off atmospheric nuclei. 
While the vast majority of these cosmic rays lose their energy through one or more hard 
QCD interactions, a small fraction can undergo new-physics interactions. With energies 
ranging up to 10 11 GeV (in the earth's frame), the incident cosmic rays are capable of 
generating center-of-mass energies above 10 5 GeV when scattering off a proton. And with 
the exception of the very highest energies, the incident cosmic ray energy spectrum is 
well-measured. 

While these weak interactions in the upper atmosphere are useless for studying short- 
lived SUSY states, such states will be produced nonetheless. And if R-parity is conserved, 
those SUSY states will eventually decay down to the LSP. Thus LSPs are produced in 
our atmosphere in standard particle physics processes, independent of cosmological or 
astrophysical assumptions 

2.1 Spectrum and Composition 

The all-particle spectrum for cosmic rays is well-measured up to energies of about 10 11 GeV 
At energies up to about 10 6 GeV the spectrum follows a simple E~ 2J power law. Around 
10 6 GeV, the so-called "knee" is hit at which point the spectrum begins falling more 
dramatically, as E~ 3 ' 2 . Finally, at energies around 10 9 GeV (the "ankle") the spectrum 
flattens out slightly to E~ 3 || |HJ. 
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The composition of the incident cosmic rays is not as well known. The primary 
spectrum is composed of a number of elements including protons, helium, iron, etc. 
For energies below the "knee" (10 6 GeV per nucleus), protons are the most abundant 
constituents, representing roughly 50% to 80% of the spectrum. At the knee, iron nuclei 
begin to dominate the spectrum; however, only qualitative details of the composition can 
be inferred using, for example, the depth of shower maximum. The data seems to indicate 
that there is a relative rise in heavier elements and then a gradual decrease again as the 
"ankle" is approached ||. Above 10 9 GeV per nucleus, lighter elements appear to be more 
abundant again, consistent with fragmentation of the heavier elements by the CMBR. 
Extractions from the data of the average particle number as a function of primary energy 
vary from 5 to 15 for different experiments ||. Given this level of ignorance, we assume 
a fixed fraction of protons (50%) at energies below 10 6 GeV, exponentially decreasing to 
a constant 1% at energies above 10 8 GeV. 

We will also assume that the remainder of the primary composition (the part that 
is not protons) has on average A nucleons/incident nucleus; we will choose A = 13 as 
a typical experimental value ||. We will model these nuclei as a collection of A loosely 
bound nucleons, each carrying E/A of the incident energy E of the entire nucleus; in 
the center-of-mass frame, the total energy of the collision, a/s, is therefore reduced by a 
factor A~ x l 2 . To be conservative, we assume that any one nuclei can only participate in 
a single hard interaction and we do not consider the interactions of the daughter nuclei, 
if any are formed, at all. We find that nuclei heavier than hydrogen are responsible for 
only about 10% of the SUSY interactions and therefore our result is almost independent 
of how we model their interactions. We believe the above choices to be very conservative. 
Looser constraints with more protons in the spectrum could allow our calculated rates 
to more than double. 



2.2 Atmospheric Cascades 

Protons [] incident on our atmosphere will usually lose their energy either by QCD or 
QED processes. Typical QED energy loss rates through brehmstrahlung and ionization 
are about 2 MeV/(g/cm 2 ) |H| so that charged particles lose only a few GeV traversing 



15 km of atmosphere. This is negligible compared to losses in hard QCD interactions. 

Energetic primary protons have a nuclear interaction length of 80g/cm 2 , correspond- 
ing to 12 QCD interactions over the depth of the atmosphere [|l(J. For simplicity, we 
only consider the first such interaction and neglect secondaries produced in the resulting 
air showers. That is, we assume that once an incident proton has had a single hard QCD 
interaction, its energy is degraded beyond the point at which it is kinematically possible 
to create SUSY particles. In effect, we are discarding the possibility that a secondary 
could participate in a SUSY interaction, a choice which again results in a conservative 
bound on the number of SUSY states produced. A complete analysis, unneccesary here, 
involves solving a set of coupled cascade equations with appropriate boundary condi- 



1 This discussion is easily generalized to neutrons or heavy nuclei. 
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tions, taking into account energy loss processes and energy-dependent cross sections. 
Our method reduces to analyzing a fixed-target pp collision with one of the protons at 
very high energies. 



3 The Calculation 

Now we present the details of our calculation. After passing through a distance dx of 
the atmosphere with local density p a tm, the flux of protons with energies between E and 
E + 5E which have undergone hard QCD interactions and are therefore "lost" for SUSY 
interactions is simply given by 

1 d ^ E K- %(E) ^ (E) 



n atm (x) dx 

where $ p is the flux of protons measured in GeV -1 cm* 2 s _1 sr -1 , n atm is the number den- 
sity of protons in the atmosphere at a depth x, and E is the energy of the incident proton 
in the earth's frame. At the energies under consideration here, the QCD pp cross-section, 
(Tp P , is roughly 100 mb, independent of energy. Over that same slice of atmosphere, the 
total number of SUSY interactions (integrated over all energies), corresponding to the 
flux of LSPs, is 

_id_^y = r dE%{E)asvsY{Eh 

n s,tm{ x ) dx J 

Here osusy represents the inclusive cross-section for pp — > Xsusy, independent of whether 
X is squarks, sleptons, etc. (osusy is a function of s = 2m p E.) This system of equations 
is simple to solve. For a spectrum of only protons, the total number of SUSY interactions 
in cm -2 s _1 sr _1 is then: 



SUSY 




where $ p and ogusY are implicitly functions of E. Accounting for the observation that 
only a fraction f p of the incident cosmic rays are protons, and using our stated assumption 
that the remaining cosmic rays contain on average A nucleons, then the above discussion 
generalizes to: 

$ SUSY = JdE (^j [f p a SVSY (E) + (1 - f p ) asusY (f ) 

where d^i/dE is the all-particle incident flux of cosmic rays. 

There is another, simpler way to understand the above calculation. Since we only 
allow each incident cosmic ray one hard interaction (and only 1 LSP per SUSY inter- 
action), we can simply calculate the probability that that one interaction will be either 
QCD or SUSY. Since <t S usy <^ c PP , the probability of a SUSY interaction is then roughly 
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ususy / 'ffpp- That probability is then integrated over the entire flux of incident cosmic 
rays to find the flux of LSPs produced. 

The only unknown quantity that remains is osusy- Since this represents a total cross- 
section, summed over all possible SUSY final states, there is considerable room for model- 
dependence. However, we can make several simplifying observations/assumptions. First, 
because the cosmic rays and the atmospheric nuclei are baryonic, strong interactions 
should dominate the SUSY production processes. Second, in most "realistic" SUSY 
models, there is an approximate degeneracy among the strongly-interacting sparticles, 
that is, among the squarks and the gluinos. These two statements allow us to greatly 
simplify the calculation by only considering superQCD interactions and by assuming that 
all squarks and gluinos have a common mass Msusy- 

Independent of the details, the parton-level cross-section for SUSY production obvi- 
ously goes like a 2 s /s. To be slightly more realistic, we will use actual calculated cross- 
sections for the range of superQCD processes. The processes that we consider are: 

(a) qq -> qq_ (c) qq -> (e) gg -> gg 

(b) gg -> qq (d) qq -> gg (f) qg -> qg 



The relevant cross-sections have been calculated and tabulated in Ref . Ill; we confine our 



calculation to tree level. To go to the proton-proton cross-section we use the CTEQ5M 



parton distribution functions [|T^] . The total pp — > SUSY cross-section is shown in Fig- 
ure |as a function of y/s for several choices of Msusy- Over the energy range of interest, 
the SUSY cross-section is dominated by qq — > qq, qq and qg — > qg. There is a subtlety 
associated with exact versus approximate degeneracy of the squarks and gluinos and so 
we consider both the case in which rriq — m- g ~ and rriq — rrig ~ Msusy- We find a 
difference of about 50% in our calculated number densities as we vary over this range, 
and so we show a range of limits. 

Once a squark or gluino has been produced, it will decay in a cascade down to the 
the LSP. Presumably, one LSP of positive charge will be produced for each of negative 
charge. Searches for the latter are more difficult and we will concentrate only on the 
former. Further, we will assume conservatively that only one positively charged LSP 
is produced per SUSY interaction, though the number can be significantly higher as 
the cascades of decays progress. One expects that these positively charged LSPs form 
superheavy hydrogen by attracting a nearby electron and that this superheavy hydrogen 
eventually (over the lifetime of the earth) bonds into a superheavy water molecule in the 
earth's oceans. 

So as a final step we calculate the concentration of superheavy water in the oceans. 
Considering the age of the oceans, to, to be roughly 3 billions years, and assuming 
that the flux of cosmic rays has remained essentially unchanged over that time period, 
the number of superheavy water molecules per usual water molecule in the oceans (the 
"anomalous concentration" , x) is 

7r m-13 ( *o \ /2.6km\ $ SU sy 

X = 7.6 x 10 — — — x — - 

\6Lrjr J \ d J cm z s 1 sr 1 
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Figure 1: Total cross-section for pp — ► SUSY as a function of y/s for several values of Msusy- 



where do is the average depth of the oceans. 

The constraints on \ come from searches for superheavy molecules in large samples 
of water. Experiments then place bounds on x as a function of the mass of the stable, 
charged particle. In our approach, this mass in unknown, though it is bounded from 
above by Msusy- The strongest bound comes the experiment of Smith et al. H, who 
find limits ranging from \ < 3.5 x 10 -30 for m^sp = 100 GeV to \ < 8 x 10 -28 for 
m LSP = 1.2 TeV. Because their bound monotonically weakens as vn^sp increases, we can 
place a very conservative bound on Mgusy by setting rriLSP — ^susy- For that case, we 
find Msusy > 230 GeV if the LSP is stable and charged. (Our bound has an uncertainty of 
roughly ±10 GeV due to a ~ 50% uncertainty that comes from the subtleties for defining 
Msusy discussed above.) These results are summarized in Figure ^| where we have shown 
the predicted anomalous concentration as a function of Msusy and the experimental 
bound of Ref. §. 



4 Conclusions 

We have shown that there exist bounds from cosmic ray production in the upper atmo- 
sphere on charged stable relics (a charged SUSY LSP in particular) which are independent 
of cosmological constraints. Under the assumption that the incident cosmic ray flux has 
remained constant over the last 3 billion years, we have calculated a conservative lower 
bound on the scale of new physics (Msusy), using nothing more than standard particle 
physics. If SUSY has a stable and charged LSP, then we can place a lower bound on the 
mass scale of the squarks and gluinos at 230 GeV. 

This procedure can easily be extended to other models of weak-scale physics. In 
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Figure 2: Anomalous concentration, of superheavy hydrogen in sea water as a function of 
the SUSY mass scale (shaded region). The bound of Ref. [|| is shown as a dashed line. The 
limits shown here assume m^sp = Msusy- The width of the shaded region represents a range 
of ±50% in the cross-section calculation due to model-dependence. 
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such a case, limits similar to those found here could be placed on the masses of new 
strongly-interacting particles. 

Importantly, these bounds will not change if the basic paradigms of cosmology at 
and below the weak scale are questioned, such as happens in models with large extra 
dimensions. For example, it is unclear how standard cosmological bounds can be used 
to constrain models such as that of Ref. [0] which predicts a light, stable top squark but 
becomes non-perturbative and higher-dimensional at the TeV scale. The bound presented 
here should hold even in these highly non-standard cases. 

Finally, if a charged particle is discovered which is stable on collider timescales, ruling 
out or verifying that it is stable on cosmological timescales will require that searches for 
superheavy water be examined again, though with larger initial samples. Unfortunately, 
the steeply falling cosmic ray spectrum requires us to go to exponentially larger sam- 
ples in order to significantly increase our sensitivity. For example, the method used in 
Ref. would require an initial heavy water sample equal to that contained in the SNO 
experiment in order to probe squark masses up to 700 GeV. While such a large-scale 
search seems unnecessary at present, any future discovery of a new stable, charged par- 
ticle might require just such an effort. Otherwise there may be no other way to study 
the stability of that state on long timescales. 
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